Approved  for  public  release;  distribution  is  unlimited. 


Optimal  Array  Geometries  for  Wideband  DOA  Estimation 


Ulkii  Baysal  and  Randolph  L.  Moses 

Department  of  Electrical  Engineering,  The  Ohio  State  University 
2015  Neil  Avenue,  Columbus,  OH  43210  USA 

Abstract 

We  consider  the  role  of  array  geometry  on  the  Direction  of  Arrival  (DOA)  estimating  performance  of  the 
array  where  the  impinging  signal  is  wideband.  We  concentrate  on  arrays  that  have  isotropic  performance. 
An  isotropic  array  is  one  whose  Cramer  Rao  bound  (CRB)  on  the  DOA  of  a  single  source  is  uniform  for  all 
angles.  We  derive  the  necessary  and  sufficient  conditions  on  the  location  of  array  elements  so  that  the  array 
is  isotropic.  Both  planar  arrays  and  three  dimensional  arrays  arc  considered.  We  also  present  several  designs 
of  isotropic  planar  and  volume  arrays  and  give  example  geometries. 


1.  Introduction 


The  number  and  location  of  the  elements  in  an  array  strongly  affects  the  Direction-of-Arrival  (DOA)  esti¬ 
mation  performance  of  the  array  system.  There  is  a  considerable  amount  of  work  done  on  the  design  of  the 
optimal  array  (optimal  in  terms  of  cost,  space,  error  variance  or  resolution  limits  etc.).  Most  of  the  emphasis  is 
devoted  to  linear  arrays  (or  combination  of  linear  arrays)  as  they  arc  simple  to  analyze,  provide  the  maximum 
aperture  when  the  number  of  elements  fixed  and  optimal  DOA  estimation  algorithms  arc  available  for  such 
arrays  [1] — [5].  One  of  the  main  problems  with  the  linear  arrays  is  the  nonuniformity  of  the  performance:  the 
DOA  estimation  performance  degrades  considerably  near  endfire.  In  this  paper,  we  concentrate  on  arrays 
that  have  uniform  performance  over  the  whole  field  of  view. 

Several  different  performance  and  design  criteria  have  been  introduced  to  be  used  in  obtaining  optimal  arrays. 
Performance  comparisons  of  some  common  array  geometries  are  presented  in  [6], [7],  In  [8],  the  authors 
introduce  a  measure  of  similarity  between  array  response  vectors  and  show  that  the  similarity  measure  can 
be  tightly  bounded  below.  The  array  with  the  highest  bound  is  optimum  in  the  sense  that  it  has  the  best 
ambiguity  resolution.  In  [9],  a  sensor  locator  polynomial  is  introduced  for  array  design.  A  polynomial  is 
constructed  using  prespecified  performance  levels,  such  as  detection-resolution  thresholds  and  Cramer-Rao 
Bounds  (CRBs)  on  error  variance,  and  its  roots  arc  the  sensor  locations  of  the  desired  linear  or  planar  array. 
In  [10],  differential  geometry  is  used  to  characterize  the  array  manifold  and  an  array  design  framework  based 
on  these  parameters  is  proposed. 
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Studies  regarding  to  arrays  that  have  uniform  performance  with  respect  to  a  certain  criterion  can  also  be 
found  in  the  literature.  In  [11],  the  asymptotic  mean  square  angular  error  is  used  to  define  an  isotropic 
array.  The  authors  derive  the  angle  CRB  for  a  single  far-held  source  and  then  derive  conditions  on  the  sensor 
locations  to  ensure  the  azimuth  and  elevation  errors  are  uncoupled  from  each  other  in  the  bound.  In  [15] 
and  [18]  the  authors  consider  conditions  on  the  array  geometry  for  which  the  single  source  azimuth  and 
elevation  CRBs  arc  uncoupled,  and  in  [15]  a  condition  which  gives  isotropic  array  performance  is  derived. 
In  [12],  some  popular  array  geometries  arc  compared  using  the  ratio  of  the  number  of  sensor  elements  to 
the  usable  aperture  as  a  performance  measure;  the  Y-shaped  array  shown  to  achieve  approximately  uniform 
angle-of-arrival  estimation  performance  with  this  criterion. 

In  this  paper,  we  study  planar  and  volume  array  geometries  that  have  isotropic  DOA  estimation  performance. 
For  planar  arrays,  the  arrays  arc  isotropic  in  the  sense  that  the  CRB  on  the  DOA  estimation  of  a  single  source 
is  uniform  for  all  source  arrival  angles  from  0  to  2tt.  For  volume  arrays  we  use  the  bound  on  the  Mean 
Square  Angular  Error  (MSAE)  as  the  criterion.  The  MS  AE  is  a  scalar  measure  of  the  error  between  true  and 
estimated  unit  bearing  vectors  pointing  towards  the  source,  and  its  bound  is  computed  from  the  CRB.  The 
array  is  said  to  be  isotropic  if  the  bound  on  the  MSAE  is  constant  for  all  azimuth  and  elevation  angles  in  [0 
2tt|  x  [-|  |],  Since  the  CRB  and  bound  on  the  MSAE  are  independent  of  any  particular  estimator  and  ML 
estimators  asymptotically  achieve  these  bounds,  they  arc  useful  criteria  for  array  design.  Our  results  apply 
to  both  narrowband  and  wideband  signals. 

The  organization  of  the  paper  is  as  follows.  In  Section  II,  we  describe  the  system  model  and  state  our 
assumptions.  Section  III  discusses  the  planar  array  scenario.  First,  we  define  the  performance  criterion  and 
isotropy  condition  for  planar  arrays,  and  later  we  state  the  necessary  and  sufficient  conditions  on  the  array 
geometry  for  isotropy.  We  also  propose  different  array  design  methods.  In  Section  IV,  we  study  the  three 
dimensional  arrays.  We  introduce  the  isotropy  condition,  and  give  the  necessary  and  sufficient  conditions  on 
the  three  dimensional  array  for  isotropic  performance.  Section  V  concludes  the  paper. 


2,  System  Model 


We  assume  an  array  of  N  identical  sensors  located  in  the  space  at  locations  r,  for  i  €  [1 .  Ar] .  We  will  consider 
both  planar  arrays  in  which  r,  =  [rXi,ryi]T  and  volume  arrays  in  which  r,  =  [rXi,  rVi,  rZi]T. 

Following  [14],  we  adopt  a  system  model  describing  a  source  impinging  on  the  array.  A  single  Tar-field 
source  s(t),  which  is  in  general  wideband,  impinges  on  the  array  from  direction  6  =  [(/>,  ip\,  4>  denoting  the 
azimuth  angle  measured  counterclockwise  from  the  x-axis  on  the  x-y  plane,  and  A  denoting  the  elevation 
angle  measured  counterclockwise  from  the  x-axis  on  the  x-z  plane.  The  noise  at  the  sensors  are  independent, 
zero  mean  Gaussian  noise,  and  independent  of  the  source  signal.  The  observation  time  T  is  partitioned  into 
K  intervals  of  length  T,i  and  a  J-point  discrete  Fourier  transform  is  applied  to  each  interval.  Assuming  Td 
is  long  enough,  we  say  the  discrete  Fourier  coefficients  arc  uncorrelated.  Then, 

xk{uj)  =  +  nk(uj),  j  =  k  =  l,...,K  (1) 

where  xk(uj),  ?ik(uj)  are  N  x  1  vectors  and  sk(u>j)  is  a  scalar.  The  elements  of  xk{t Oj),  nk(iOj)  and  sk(ajj) 

are  the  discrete  Fourier  coefficients  of  the  sum  of  the  sensor  outputs,  the  noise  and  the  signal  source  at  the 


discrete  frequency  u3,  respectively.  Ag(uj)  is  given  by 


Agiujj)  =  [e^d^d\ejuJ^e\  . . . ,  eM<bv(0)]T  (2) 

where  dk(0)  =  u  erfc  is  the  propagation  delay  associated  with  the  kth  sensor,  c  is  the  speed  of  propagation 
and  u  denotes  the  unit  bearing  vector  pointing  towards  the  signal  source.  For  a  planar  signal  arriving  from 
angle  (p 

u  =  [cos(0),  sin(c/>)]T  (3) 

whereas  for  a  three-dimensional  case  where  a  signal  arrives  from  azimuth  angle  <f>  and  elevation  angle  ©, 

u  =  [cos(0)  cos(V’),  sin(</>)  cos^ip),  sin(V’)]T  (4) 


3.  Planar  Arrays 


In  this  section  we  consider  the  special  case  of  a  planar  array  with  elements  at  locations  r,  =  [rXi ,  ryi]T.  The 
array  is  used  to  estimate  the  DOA  of  a  wideband  signal  s(t)  which  is  coplanar  with  the  array.  Here  the  signal 
arrives  at  an  angle  (t>. 


3. 1  Single-Source  CRB 


For  the  system  model  described  in  Section  II  and  under  the  planar  array-coplanar  signal  assumption,  the 
CRB  for  the  source  DOA  estimate  is  given  by  [14]: 


CRB((j)) 


(5) 


Here,  A $  =  ^4,  P\  =  I  —  A(An Ay  lArf  is  the  projection  matrix  onto  the  subspace  orthogonal  to  the 
column  space  of  A,  Rs(u:j)  is  the  cross-spectral  density  matrix  of  the  impinging  signals  at  frequency  t Oj  and 
©  denotes  the  Hadamard  product. 


For  the  model  given  by  (1)  and  (2),  we  can  write  the  CRB  as 


CRB{4>) 
G(B,  (/>) 
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G(B,4 i) 


-t  /  2ii  A 
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Pj  1  - 
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cos(0) 
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(6) 

(V) 


(8) 


where  rc  is  the  centroid  of  the  array,  i.e., 


(9) 

i=  1 

Pj  is  the  signal  power  and  rij  is  the  noise  power  at  frequency  interval  j. 

We  see  that  the  CRB  is  a  product  of  a  term  G(B ,  0)  that  depends  only  on  the  source  DOA  and  the  array 
geo  me  try  and  a  term  that  depends  on  source  and  noise  powers  as  a  function  of  frequency.  This  is  an  important 
property,  because  the  impact  of  the  array  geometry  on  the  CRB  is  the  same  regardless  of  whether  the  source 
spectrum  in  narrowband  or  broadband,  and  regardless  of  the  source  signal  and  noise  spectral  densities.  Thus, 
the  results  that  follow  apply  to  a  broad  class  of  array  signal  processing  scenarios.  Moreover,  the  CRB  depends 
on  the  array  geometry  only  through  the  matrix  B,  which  is  the  2x2  “covariance”  of  the  array  points.  Thus,  any 
two  geometries  that  have  the  same  covariance  matrix  B  will  have  identical  single-source  CRB  performance. 

We  note  that  the  array  performance  criterion  we  have  chosen  does  not  take  into  account  potential  array 
ambiguities  that  arise  when  the  array  manifold  from  two  different  DOAs  arc  close  to  one  another  (see  [8] 
for  a  discussion  of  this  topic).  For  wideband  arrays,  ambiguities  arc  not  much  of  a  problem  because  the 
frequency  diversity  eliminates  most  or  all  DOA  ambiguities. 


3.2  Geometric  Interpretation 


It  is  possible  to  write  G{B ,  0)  explicitly  in  the  following  form. 


N 


G(B ,  0)  =  ^  1 1 r- * 1 1 2  sin2  (  0  —  tan(  — 

"  "  V  V  rxi 


i=  1 


l 
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r  N 

Ei 

i= 1 


sin 


(10) 


The  following  geometrical  interpretation  can  be  derived  from  the  above  expression  for  the  geometry- 
dependent  term  G(B ,  0).  Project  the  N  sensor  points  onto  a  line  that  is  perpendicular  to  the  line  passing 
through  the  origin  and  orthogonal  to  the  direction  vector  u.  (See  Figure  1).  Then  G(B.  0)  is  the  sample 
variance  of  these  projected  points. 

Note  also  that  G(B.  0)  can  be  interpreted  as  a  measure  of  the  beamwidth  of  the  array.  Consider  a  fixed 
frequency  cco  and  corresponding  wavelength  Ao  =  If  we  choose  the  delay-and-sum  array  weights  to 
steer  the  beam  at  an  angle  0q.  then  resulting  complex-valued  array  response  at  angle  0  at  this  frequency  is 
given  by 


w4>0{(t))  =  aH(0)a{0  o)  (11) 

a{0)  =  [ej^uT-ri,...,ej%uT-rN]T  (12) 

The  array  gain  WA0(ci)|  can  be  approximated  by  a  quadratic  about  the  point  0  =  0q.  A  second  order  Taylor 
series  expansion  of  |FF</,0(^>)|  yields 

\Wto(0)\~N-^(^J  G(B,  0o)(0  —  0o)2  (13) 


Figure  1:  Geometrical  interpretation  of  G(B,d>)  as  the  sample  variance  of  the  array  locations  projected  onto 
the  line  orthogonal  to  u 


Using  (13),  the  half-power  beamwidth  of  the  array  is  given  by  fadB  =  <fto  +  ^  q{B  <j>0)  •  Thus,  \JG(  B.  cyi) 
is  inversely  proportional  to  the  half-power  beamwidth  of  the  array. 


3.3  Isotropic  Planar  Arrays 

From  equation  (8)  we  see  that  G(B,  <fr),  and  also  the  single-source  CRB,  is  independent  of  translation  of  the 
array  element  locations,  so  without  loss  of  generality  we  can  assume  that  the  array  is  centered  at  the  origin, 
or  equivalently,  rc  =  [0,  0]T.  Under  this  assumption,  B  in  (8)  simplifies  to 

1  N 

B=N^2riri  (14) 

i= 1 

We  are  interested  in  planar  geometries  whose  single-source  CRB  is  independent  of  signal  arrival  angle  <fi. 
We  refer  to  such  arrays  as  isotropic  arrays.  The  following  result  characterizes  the  set  of  all  isotropic  arrays: 

Theorem  1: 

(a)  An  N  element  planar  array  which  is  centered  at  the  origin  (rc  =  [0,  0]T)  and  represented  by  the  2-by-2 
matrix  B  =  ^  YliLi  rirf-  where  r*  =  [rXi,ryiY  is  the  location  of  the  ith  sensor,  is  isotropic  if  and  only  if 

B  =  kl2  (15) 


where  k  is  any  positive  constant  and  1-2  denotes  the  2x2  identity  matrix. 

(ft)  If  we  parameterize  the  array  geometry  by  an  N  x  1  complex-valued  vector  g  given  by  g  =  [V/i , . . . ,  g^]T 
where 

9i  =  rXi+jryi  (16) 

then  condition  stated  in  (15)  is  equivalent  to 

gTg  =  o  (17) 


(c)  If  the  array  is  isotropic,  then 


G(B,</>)  =  k 


1=1 


(18) 


Proof:  See  [19], 


3.4  Planar  Array  Design  Examples 

In  this  section  we  present  four  design  methods  for  generating  planar  arrays  and  give  examples  for  each 
method. 

Circularly  Symmetric  Geometries: 

The  immediate  solution  where  both  equations  (14)  and  (15)  arc  satisfied  is  circularly  symmetric  geometries 
or  any  superposition  of  circularly  symmetric  geometries,  where  we  define  a  circularly  symmetric  geometry 
as  one  in  which  N  >  3  sensors  arc  equally-spaced  on  a  circle  that  has  a  radius  greater  than  0.  A  single  sensor 
located  at  the  origin  is  also  in  this  class.  An  example  geometry,  where  a  4  element  and  a  3  element  circularly 
symmetric  geometry  arc  combined  to  form  a  7  element  array,  is  given  in  Figure  2.  Notice  that  radii  of  the 
subarrays  may  be  equal  or  unequal. 


Figure  2:  An  isotropic  planar  array  obtained  as  the  superposition  of  a  4  element  and  a  3  element  circularly 
symmetric  array. 


Rotated  Geometries: 

A  second  class  of  isotropic  planar  arrays  can  be  obtained  as  follows.  Consider  an  element  subarray,  where 
N  is  even,  with  arbitrary  sensor  locations.  Define  the  origin  as  the  centroid  of  these  points,  and  define  a 
second  ^  element  subarray  by  rotating  the  first  elements  by  either  90°  or  —90°.  Consider  the  iV-element 
array  formed  by  the  superposition  of  those  two  subarrays.  Note  that  the  centroid  of  the  rotated  elements 
is  the  origin,  as  is  the  centroid  of  the  JV-element  array.  Moreover,  in  complex  coordinates,  the  iV-element 

array  vector  g  can  be  written  as  g  =  gTN  e3  2 ,  g1^,  ] T,  where  <7  jv  parameterizes  the  element  array.  Then 

~2  T  2 

gT g  =  g^gN  +  e311  grN  g<v  =  0,  hence  any  geometry  formed  following  this  procedure  is  isotropic.  An 

T  2  "2"  2 

example  array,  generated  by  randomly  selecting  the  locations  of  the  first  three  elements,  is  shown  in  Figure  3. 


90 

1 


Figure  3 :  A  6  element  isotropic  array  (circles  and  stars)  formed  from  rotating  a  randomly  generated  3-element 
subarray  (stars)  by  90°. 


More  generally,  one  can  take  any  arbitrary  subarray  of  size  N/m  and  combine  its  0,  2n/m, . . . ,  2n  (m—  l)/m 
rotated  versions  about  any  arbitrary  point  (not  just  the  center  of  gravity);  the  resulting  N  elements  form  an 
isotropic  array.  This  array  is  also  a  superposition  of  m  N/m-e lenient  circularly  symmetric  geometries, 
though,  and  covered  in  the  previous  example. 

Completion  of  Arbitrary  (N-2)-element  Arrays: 

It  is  possible  to  obtain  an  N  element  isotropic  array  from  an  arbitrary  (N  —  2)  element  array  by  adding 
two  elements.  Assume  that  locations  of  the  first  (N  —  2)  sensors  arc  given  and  let  rjv_  1  and  r/y  denote  the 
locations  of  the  remaining  two  sensors.  Setting  9i  =  0  anc*  9T 9  =  0  (see  (16))  gives 

N- 2 

9N- 1  +  9N  =  —  ^2  9l 
1=1 

N—2 

2  1  2  \  "  2 

9n—  1  +  9n  ~  ~  Si 
1=1 

The  solutions  to  (19)  and  (20)  uniquely  determine  the  locations  of  the  last  two  sensors  so  that  the  resulting 
A' -clement  array  is  isotropic.  An  example  geometry  formed  with  this  method  is  shown  in  Figure  4. 


(19) 

(20) 


90 
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Figure  4:  A  7-element  isotropic  array  formed  from  an  arbitrary  5 -element  subarray  (stars)  by  adding  2 
elements  (circles). 


X-shaped  Isotropic  Arrays: 

We  can  combine  two  X-shaped  geometries  so  that  the  resulting  8-element  array  is  isotropic.  An  X-shaped 
geometry  is  a  set  of  four  sensors  with  radii  ||  r||  and  angles  ±a,  ±7r  —  a.  It  can  be  shown  that  any  pair  of 
X-shaped  geometries  with  parameters  (||rj||,  aii),  i  =  1 ,  2  that  satisfy  ||n ||2  cos(2cri)  +  ||r2 1|2  cos(2a2)  =  0 
is  an  isotropic  array. 

A  special  case  is  the  superposition  of  X-shaped  geometries  arc  those  whose  elements  lie  along  parallel  two 
lines.  For  two  superimposed  X-shaped  geometries,  we  constrain  ||ri||  cos(cr)  =  |J ?™2 1|  cos(cr).  For  example, 
the  eight-element  isotropic  array  with  elements  having  ./;- val  ucs  of  ±1  is  shown  in  Figure  5. 


90 
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Figure  5:  An  8-element  isotropic  array  formed  by  combining  two  X-shaped  geometries. 


4.  Three  Dimensional  Arrays 


In  this  section,  we  consider  an  array  that  has  elements  located  in  R3  and  is  used  to  estimate  the  DOA  of  a 
single  wideband  far-held  source.  The  source  direction  is  parameterized  by  9  =  [< j),ip}T,  where  6  €  [0,  2i r) 
and  f  G  [— f ,  f  ]  denote,  respectively,  the  azimuth  and  the  elevation  of  the  source.  The  single  source  CRB 
for  this  scenario  is  a  2-by-2  matrix, 


CRB(9) 
G(B,e ) 


Jg(u)TBJg(u ) 


3  =  1 


, 

— Pj  I  1  - 
n- 


71  n 


pjN  +  rij 


-l 


(21) 

(22) 


where  B  is  given  by  (8)  and  Jq(u)  is  the  3-by-2  Jacobian  matrix  of  the  u  given  in  (4),  and  is  equal  to 


—  sin (<fi)  cos(^) 

—  cos((j>)  sin(r/>) 

du 

dip 

du 

dip 

= 

cos (f)  cos(r/>) 

—  sin(</>)  cos(^) 

0 

cos(V’) 

(23) 


Once  again,  wee  see  that  the  CRB  is  the  product  of  a  term  that  depends  only  on  the  array  geometry  of  source 
angle,  and  a  term  that  depends  only  on  the  signal  and  noise  parameters.  We  also  see  that  the  CRB  depends 
on  array  geometry  only  through  the  3x3  array  covariance  matrix  B. 


4.1  Performance  Criterion 

Estimating  the  azimuth  and  elevation  corresponding  to  the  DOA  of  the  signal  is  equivalent  to  estimating  the 
vector  u.  A  direction  6  =  [©,  f>]T  uniquely  specifies  u  via  (4)  and  an  estimate  9  =  [f,  C]  uniquely  specifies 
a  vector  u.  Let  S  be  the  angle  between  the  vectors  u  and  u;  since  both  u  and  u  are  unit  vectors,  we  have 

cos  (5)  =  uT  ■  u  (24) 


The  Mean-Square  Angular  Error  (MSAE)  is  introduced  in  [16]  as  a  scalar  measure  of  estimator  performance 
in  estimating  a  geometrical  vector  and  defined  as  the  expectation  of  52.  The  MSAE  enjoys  two  desired 
properties:  it  is  independent  of  the  choice  of  the  reference  coordinate  frame  and  it  does  not  suffer  from 
the  singularity  inherent  in  spherical  coordinates  as  ip  — >  +|.  The  lower  bound  of  the  MSAE  provides  a 
performance  criterion  for  a  set  of  estimators  that  satisfies  certain  mild  conditions,  similar  to  those  needed  for 
the  CRB.  A  derivation  for  the  lower  bound  of  the  MSAE  and  a  detailed  discussion  of  the  conditions  for  the 
applicability  and  tightness  of  the  bound  can  be  found  in  [17].  Assuming  M  is  the  number  of  observations, 
the  asymptotic  normalized  MSAE  is  defined  as 

MSAE oo  =  lirn  MES2  (25) 

M—>  oo 

and  when  spherical  coordinates  arc  used,  (25)  is  bounded  below  as 

MSAEoo  >  MSAEb  =  cos 2('tJj)CRB00(f)  +  CRB 0O(V’)  (26) 

where  CRB ^  is  the  asymptotic  CRB.  For  a  geometrical  interpretation  of  the  MSAE#  see  [11]. 


4. 2  Isotropic  Three  Dimensional  Arrays 


We  adopt  MSAE#  as  a  performance  criterion,  and  define  a  three  dimensional  array  to  be  isotropic  if  the 
associated  MSAE#  is  constant  for  all  [</>,  f>]T  €  [0,  27t)  x  [— |],  The  following  theorem  defines  the  set  of 
all  isotropic  three  dimensional  arrays: 

Theorem  2: 

(a)  An  N  element  array  which  is  centered  at  the  origin  (rc  =  [0,  0,  0]  7  )  and  represented  by  the  3-by-3  matrix 
B  =  Eil l  nrj,  where  n  =  [rXi,  ryt ,  r2J  7  is  the  location  of  the  ith  sensor,  is  isotropic  if  and  only  if 

B  =  kh  (27) 


where  k  is  any  positive  constant  and  /  >  is  the  3-by-3  identity  matrix. 

(b)  If  (27)  holds 

G(B,(t>)  =  kI2  (28) 


Proof:  See  [19], 

We  remark  that  (27)  is  also  necessary  and  sufficient  for  the  2x2  CRB  matrix  to  be  independent  of  source 
arrival  angle.  In  this  case,  the  CRB  in  (22)  is  diagonal,  and  if  it  is  scaled  to  remove  the  latitudinal  scaling 
of  azimuth,  the  CRB  in  the  azimuthal  and  elevation  directions  are  equal.  That  is,  an  uncertainty  ellipse  in 
spherical  angle  is  a  circle  whose  radius  is  independent  of  source  arrival  angle. 

In  [11],  the  authors  give  sufficient  conditions  on  the  array  geometry  so  that  MSAE/j  is  independent  of 
the  source  signal  DOA.  The  above  theorem  extends  their  results  by  proving  that  these  conditions  are  also 
necessary. 


4.3  Three  Dimensional  Isotropic  Array  Example 

Analogously  to  the  planar  case,  it  can  be  shown  that  arrays  formed  by  placing  the  sensor  elements  at  vertices 
of  any  regular  polyhedron1,  or  a  superposition  of  such  arrays,  result  in  three  dimensional  isotropic  arrays. 
We  conjecture  that  the  result  also  holds  for  the  13  semi  regular  polyhedra. 

It  can  also  be  shown  that  three  elements  can  be  added  to  an  arbitrary  (N  —  3)-element  subarray  to  make  the 
resulting  iV-element  array  isotropic.  Simple  design  formulas  for  the  locations  of  the  three  elements,  similar 
to  equations  (19)-(20),  can  be  derived. 


5.  Conclusion 


In  this  paper,  we  studied  planar  and  three  dimensional  arrays  that  have  isotropic  performance.  For  planar 
arrays,  we  adopted  the  single  source  wideband  Cramer-Rao  bound  as  the  performance  criterion  and  derived 


1  there  are  five  regular  polyhedra:  the  tetrahedron,  cube,  octahedron,  dodecahedron,  and  icosahedron. 


the  necessary  and  sufficient  conditions  on  the  location  of  sensor  elements  so  that  the  CRB  is  constant  for  all 
arrival  angles.  These  conditions  are  valid  regardless  of  the  source’s  frequency  spectrum.  We  presented  four 
methods  to  design  isotropic  planar  arrays. 

For  three  dimensional  arrays,  we  chose  the  asymptotic  Mean  Square  Angular  Error  as  a  measure  for  array 
isotropy.  We  derived  necessary  and  sufficient  conditions  on  the  array  geometry  that  ensure  that  the  MSAE/j 
is  independent  of  source  azimuth  and  elevation  arrival  angle.  When  these  conditions  are  satisfied,  the  azimuth 
and  elevation  are  uncoupled  in  the  CRB,  and  the  CRB  is  independent  of  source  signal  arrival  angle. 

When  designing  isotropic  arrays,  an  important  practical  issue  that  should  be  taken  into  account  is  the  minimum 
allowable  distance  between  sensors.  We  have  assumed  the  noise  components  at  the  sensors  arc  independent 
of  each  other;  this  assumption  is  violated  when  the  sensor  distances  become  small  (see  [3]).  Our  design 
methods  do  not  guarantee  that  the  resulting  sensors  are  sufficiently  well-separated,  so  if  the  designs  produce 
closely-spaced  sensors,  they  should  be  modified  accordingly. 
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